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A search for the decays B® — > )i + fx~ and B° — > fi + fi~ is performed with about 37pb _1 of 
pp collisions at yfs — 7TeV collected by the LHCb experiment at the Large Hadron Collider 
at CERN. The observed numbers of events are consistent with the background expectations. 
The resulting upper limits on the branching fractions are B(B® — > /j, + fi~) < 5.6 x 10 -8 and 
B{B° -> fi + fi~) < 1.5 x 10" 8 at 95% confidence level. 

1 Introduction 

Within the Standard Model (SM) the B® — > and B® — > flavour changing neutral 

current transitions are rare as they occur only via loop diagrams and are helicity suppressed Q 
Their branching fractions (BF) predicted by the SM 1 are B{B° S -> fJ, + fT) = (0.32 ±0.02) x 10~ 8 
and B(B° ->■ jU+aO = (0-010 ± 0.001) x 10~ 8 . 

In New Physics models, their BF can be significantly enhanced, although in some models 
they are lowered. For instance, within the Minimal Supersymmetric SM^, the BF can get 
contributions oc tan 6 f3/M\, where Ma denotes the pseudoscalar Higgs mass and tan j3 the ratio 
of Higgs vacuum expectation values. 

The best published 95% C.L. limits come from the DO collaboration 131 (6.1 fbT 1 ), B{B° S ->■ 
< 5.1 x 10~ 8 , and from the CDF collaboration 4 (2 fb -1 ), B{B° S -> fi + fi~) < 5.8 x 10" 8 
and B(B° — > /i + ^ _ ) < 1.8 x 10~ 8 . CDF has also presented preliminary results 1 ^ 1 with 3.7 fb _1 , 
that lower the limits to B(B° S -> fJ, + /J,~) < 4.3 x 10" 8 and B(B° -»■ /i + // - ) < 0.76 x 10~ 8 . 

The measurements presented in these proceedings use about 37pb _1 of integrated luminosity 
collected by LHCb between July and October 2010 at y/s = 7TeV. 

2 The LHCb detector 

The LHCb detector- 6 is a single-arm forward spectrometer with an angular coverage from about 
lOmrad to 300mrad in the bending plane. It consists of a vertex locator, a warm dipole magnet, 
a tracking system, two RICH detectors, a calorimeter system and a muon system. 

Track momenta are measured to a precision of 5p/p = 0.35 (0.5) % at 5 (100) GeV/c. The 
RICH system provides charged hadron identification in a momentum range 2-100 GeV/c. Typ- 
ically, kaon identification efficiencies of over 90% can be attained for a ir — > K fake rate below 
10%. The calorimeter system identifies high transverse energy (-Br) hadron, electron and photon 



"In this proceedings the inclusion of charge-conjugate states is implicit. 



candidates and provides information for the trigger. The muon system provides information for 
the trigger and muon identification with an efficiency of ~ 95 % for a misidentification rate of 
about 1-2 % for momenta above 10 GeV/c. 

LHCb has a two-level flexible and efficient trigger system exploiting the finite lifetime and 
large mass of heavy flavour hadrons to distinguish them from the dominant light quark processes. 
The first trigger level is implemented in hardware while the second trigger level is software 
implemented. The forward geometry allows the LHCb first level trigger to collect events with 
one or two muons with Pt values as low as 1.4GeV/c for single muon and -Pt(mi) > 0.48GeV/c 
and Pt{^2) > 0.56GeV/c for dimuon triggers. The -E/p threshold for the hadron trigger varied 
in the range 2.6 to 3.6 GeV. 

3 Analysis Strategy 

The search for £?P% — > [A + fJ,~ at LHCb is described in detail in Ref. 7 . Most of the background 
is removed by selection cuts, keeping ~ 60% of the reconstructed signal decays. Then each 
event is given a probability to be signal or background in a two-dimensional probability space 
defined by the dimuon invariant mass and a multivariate discriminant likelihood, the Geometrical 
Likelihood (GL)LL2LLH The compatibility of the observed distribution of events with a given BF 
hypothesis is computed using the CL S method^. The number of expected signal events is 
evaluated by normalizing with channels of known BF to ensure that knowledge of the absolute 
luminosity and bb production cross-section is not required. 

3.1 Event selection 

The selection consists of loose requirements on track separation from the interaction point, decay 
vertex quality and compatibility of the reconstructed origin of the B meson with the interaction 
point. The selection cuts were defined in simulation before starting data analysis. Events passing 
the selection are considered B9 s — > candidates if their invariant mass lies within 60 MeV/c 2 

of the nominal B9 \ mass. A similar selection is applied to the normalization channels, in order 
to minimize systematic errors in the ratio of efficiencies. Assuming the BF predicted by the 
SM, 0.04 (0.3) B9s —t- fj,^ pT events are expected after all selection requirements. There are 343 
(342) B?^ —t- fi> + fJ.~ candidates selected from data in the B® (B°) mass window. 

The dominant background after the B9-, — > /U + /i~ selection is expected to be bb — > nfiX® 
This is confirmed by comparing the expected yield and the kinematic distributions of the side- 
band data with a bb — > fi^X MC sample. The muon misidentification probability as a function 
of momentum obtained from data using Kg — > tt + tt^, A — > pir~ and 4> — > K + K~ decays is in 
good agreement with MC expectations. It was found that the background from misidentified 
B® d —7- h + h ~ is negligible for the amount of data used in this analysis. 

3.2 Signal and background likelihoods 

The discrimination of the signal from the background is achieved through the combination of 
two independent variables: the GL and the invariant mass. The invariant mass in the search 
regions (±60MeV/c 2 around the B9 g) masses) is divided into six equal- width bins, and the GL 
into four equal-width bins distributed between zero and one. 

The GL combines information related with the topology and kinematics of the event as the 
B9 S \ lifetime, the minimum impact parameter of the two muons, the distance of closest approach 
of the two tracks, the impact parameter and px and the isolation of the muons with respect 
to the other tracks of the event. These variables are combined using the method described in 



RefJ 10 ! 11 ! The expected GL distribution for signal events is flat, while for background events it 
falls exponentially. 

The GL variable is defined using MC simulations but calibrated with data using B9. — > 

h + h~ selected as the signal events and triggered independently on the signal in order to avoid 
trigger bias. The number of B9\ — )■ h + h~ events in each GL bin is obtained from a fit to the 
inclusive mass distribution. 

The B9s — > mass resolution is estimated from data via two methods. The first uses 

an interpolation to M B o between the measured resolutions for cc resonances (J/fy>, ip(2S)) and bb 
resonances (T(1S), T(2S), T(35)) decaying into two muons. Both methods yield consistent re- 
sults and their weighted average, a = 26.7±0.9 MeV/c 2 , is taken as the invariant mass resolution 
for both B° and B® decays. 

The prediction of the number of background events in the signal regions is obtained by 
fitting an exponential function to the \i\x mass sidebands independently in each GL bin. The 
mass sidebands are defined in the range between M B o ±600 (1200) MeV/c 2 for the lower (upper) 

two GL bins, excluding the two search regions (M B o ± 60MeV/c 2 ). 

O) 

3. 3 Normalization 

The number of expected signal events is evaluated by normalizing with three channels of known 
BF, B + — > J/ipK + , Bg —?• J/ipfj) an d B° — > K + ir~. The first two decays have similar trigger and 
muon identification efficiency to the signal but a different number of final-state particles, while 
the third channel has the same two-body topology but cannot be efficiently selected with the 
muon triggers. 



For a given normalization channel, B(B? s — > ) can be calculated as: 
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where a B o .„„ denotes the normalization factor, f B o the probability that a 6-quark fragments 

(s) W (s) 

into a B9 s and / mrm the probability that a 6-quark fragments into the 6-hadron relevant for the 
chosen normalization channel with BF B norm . The reconstruction efficiency (g REC ) includes the 
acceptance and particle identification, e SEL l REC denotes the selection efficiency on reconstructed 
events and e TRIG l SEL i s the trigger efficiency on selected events. The ratios of reconstruction 
and selection efficiencies are estimated from simulations and checked with data, while the ratios 
of trigger efficiencies on selected events are determined from data^. 

The normalization factors calculated using the three channels give compatible results; the 
final normalization factors are weighted averages and read a B o = (8.6 ± 1.1) x 10 -9 and 

flflo^ = (2-24 ± 0.16) x 10- 9 . 



4 Results 



For each of the 24 bins, the expected number of background events is computed from the fits to 
the invariant mass sidebands. For a given BF hypothesis, the expected number of signal events 
is computed using the normalization factors and the signal likelihoods. For each BF hypothesis, 
the compatibility of the expected distributions with the observed distribution is evaluated using 
the CL S method 131 . 

The observed distribution of CLg as a function of BF hypothesis can be seen in Fig. [TJ The 
expected distributions of possible values of CL S assuming the background-only hypothesis are 




Figure 1: (a) Observed (solid curve) and expected (dashed curve) CL S values as a function of B{B® — > fi~). 
The green shaded area contains the ±ler interval of possible results compatible with the expected value when 
only background is observed. The 90 % (95 %) C.L. observed value is identified by the solid (dashed) line, (b) the 

same for B(B° —> /i + /i~). 



also shown in the same figure as a green area covering the region of =_1<7 of background-only 
compatible observations. The uncertainties in the signal and background likelihoods and normal- 
ization factors are used to compute the uncertainties in the background and signal predictions. 
The upper limits read: 

B{B° S -> /.i + /0 < 4.3 (5.6) x 1(T 8 at 90% (95%) C.L., 
B(B° < 1.2 (1.5) x 1(T 8 at 90% (95%) C.L., 

while the expected values of the limits are B(B® — > /U + /i~) < 5.1 (6.5) x 10~ 8 and B(B° — > 
H + H~) < 1.4 (1.8) x 10~ 8 at 90% (95%) C.L. The limits observed are similar to the best 
published limits for the B® — > fi + decay and more restrictive for the B° — > /i + /U~ decay. 
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